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Appendices 


INTRODUCTION 


This  program  seeks  to  obtain  a  fundamental  understanding  of  the  chemical 
mechanisms  by  which  enzymes  repair  damaged  DNA,  and  to  use  this  information  to 
design  small  molecule  inhibitors  of  these  enzymes.  The  driving  force  for  these  efforts  is 
the  recognition  that  the  effectiveness  of  cancer  chemotherapy  regimes  is  intimately 
connected  to,  and  in  some  cases  directly  relies  on,  DNA  damage  repair  pathways.  A  more 
sophisticated  understanding  of  the  roles  of  DNA  damage  repair  in  the  pharmacology  of 
DNA  replication  inhibitors  will  allow  for  the  design  of  better  treatments  against  breast 
and  other  cancers. 

Studies  during  this  period  have  focused  primarily  on  Task  1  in  the  original 
Statement  of  Work,  and  this  task  has  been  essentially  completed.  We  have  also  explored 
the  question  of  how  uracil  DNA  glycosylase  (UDG)  binds  substrates  with  such  high 
affinity,  with  the  goal  of  using  this  knowledge  to  synthesize  better  rationally  designed 
inhibitors  of  the  enzyme.  Finally,  we  have  begun  work  on  Task  2  by  developing  a  new 
high  throughput  screen  for  UDG  that  will  allow  us  to  identify  inhibitors  from  large 
chemical  libraries.  This  reporting  period  has  produced  two  publications  in  peer-reviewed 
journals. 

BODY 


As  described  in  Specific  Aim  1  of  the  proposal,  we  have  now  designed, 
synthesized  and  screened  mechanism-based  inhibitors  of  UDG  (see  Appendix, 
publication  1).  This  work,  which  is  in  press  in  the  journal  Bioorganic  Chemistry,  has 
uncovered  two  useful  inhibitors  that  are  based  on  the  nature  of  the  two  transition-states 
for  the  reaction.  We  have  screened  a  directed  library  of  uracil  analogues  for  those  that 
act  most  potently  against  the  enzyme,  and  we  have  gone  on  to  tether  these  analogues  to 
the  1-azadeoxyribose  sugar,  as  described  in  the  original  proposal.  The  molecules  that 
have  emerged  from  these  studies  are  the  most  potent  inhibitors  of  UDG  yet  discovered 
(K-,  values  in  the  range  100  pM  to  14 
nM). 

We  have  moved  on  to  investigate 
what  features  of  the  DNA  substrate  are 
important  for  high  affinity  binding.  This 
work,  which  has  been  published  in  the 
journal  Biochemistry  (see  Appendix, 
publication  2),  has  revealed  that  tight 
binding  molecules  can  be  constructed  if 
the  uracil  substrate  is  presented  to  the 
enzyme  in  the  context  of  a  base  pair  that 
lacks  hydrogen  bonds. 

Finally,  we  have  developed  a 
robust  high  throughput  assay  that  allows 
the  detection  of  abasic  sites  in  DNA 
generated  by  UDG  (Fig.  1).  This  assay 


Figure  1.  High  throughput  molecular  beacon  assay 
that  detects  abasic  sites  generated  by  the  action  of 
UDG.  The  assay  relies  on  the  action  of  the  abasic 
site  endonuclease  (HAP),  which  cleaves  DNA  at 
abasic  sites.  The  fluorescent  strand  (Hex) 
dissociates  after  HAP  cleavage  resulting  in  a  20-fold 
increase  in  fluorescence  due  to  release  from  the 
quencher  (DAB).  The  assay  has  been  implemented 
in  96-well  and  384-well  formats 
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will  allow  rapid  screening  of  chemical  libraries  as  we  describe  in  Specific  Aim  2  of  the 
proposal.  This  screening  goal  will  be  performed  in  the  next  funding  period. 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  Synthesis,  screening  and  characterization  of  uracil  analogs  that  bind  to  the  UDG- 
uracil  pocket  (SOW,  Task  1). 

•  Performed  covalent  linkage  of  uracil  analogues  to  the  1-azadeoyribose  sugar 
scaffold  and  determined  the  binding  affinities  of  the  tethered  molecules  to  UDG. 

•  Elucidated  the  molecular  determinants  for  specific  recognition  of  substrates  by 
UDG  (SOW,  Task  1). 

•  Developed  a  novel  high-throughput  kinetic  assay  for  screening  inhibitors  of  UDG 
(SOW,  Task  2) 


REPORTABLE  OUTCOMES: 

This  work  has  resulted  in  two  published  manuscripts: 

1.  Jiang,  Y.  L.  Cao,  C.  ,  Stivers,  J.  T.,  Song,  F.  Ichikawa,  Y.  (2004)  The  Merits  of 
Bipartite  Transition-State  Mimics  for  Inhibition  of  Uracil  DNA  Glycosylase. 
Bioorg  Chem  In  press. 

2.  Krosky,  D.  J.,  Schwarz,  F.  P.  and  Stivers,  J.  T.  (2004)  Linear  Free  Energy 
Correlations  for  Enzymatic  Base  Flipping:  How  Do  Damaged  Base  Pairs 
Facilitate  Specific  Recognition?  Biochemistry  43, 4188-4195. 

CONCLUSIONS: 

The  long-term  goal  of  this  research  is  to  increase  the  effectiveness  of  5-FU 
chemotherapy  through  the  selective  targeting  of  BER  repair  enzymes  such  as  UDG.  We 
have  made  significant  progress  towards  this  goal  in  this  funding  period  through  the 
synthesis  of  very  potent  inhibitors  of  the  enzyme.  The  approaches  that  we  have  pioneered 
here  for  UDG  may  be  generally  applicable  to  other  DNA  glycosylases,  and  may  provide 
general  tools  for  targeting  DNA  repair  pathways.  The  next  funding  period  should  bring 
significant  progress  towards  the  development  of  cell  permeable  inhibitors  of  UDG. 

REFERENCES:  None 


APPENDICES:  Two  manuscripts  are  included  (see  above). 
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Abstract 

The  glycosidic  bond  hydrolysis  reaction  of  the  enzyme  uracil  DNA  glycosylase  (UDG)  oc¬ 
curs  by  a  two-step  mechanism  involving  complete  bond  breakage  to  the  uracil  anion  leaving 
group  in  the  first  step,  formation  of  a  discrete  glycosyl  cation-uracil  anion  intermediate,  fol¬ 
lowed  by  water  attack  in  a  second  transition-state  leading  to  the  enzyme-bound  products  of 
uracil  and  abasic  DNA.  We  have  synthesized  and  determined  the  binding  affinities  of  unimo- 
lecular  mimics  of  the  substrate  and  first  transition-state  (TS1)  in  which  the  uracil  base  is  co¬ 
valently  attached  to  the  sugar,  and  in  addition,  bimolecular  mimics  of  the  second  addition 
transition  state  (TS2)  in  which  the  base  and  sugar  are  detached.  We  find  that  the  bipartite 
mimics  of  TS2  are  superior  to  the  TS1  mimics.  These  results  indicate  that  bipartite  TS2  inhib¬ 
itors  could  be  useful  for  inhibition  of  glycosylases  that  proceed  by  stepwise  reaction  mecha¬ 
nisms. 

©  2004  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

As  a  DNA  repair  enzyme,  uracil  DNA  glycosylase  catalyzes  the  removal  of 
uracil  bases  from  DNA  that  may  arise  from  the  deamination  of  the  normal  base 
cytosine,  or  by  misincorporation  of  dUTP  into  DNA  during  DNA  replication 
[1],  Although  the  genome  protective  role  of  this  enzyme  is  well-established,  re¬ 
moval  of  uracil  is  also  important  or  essential  for  the  life  cycle  of  several  viruses, 
including  pox-,  herpes-  and  cytomegalo-,  and  thus  inhibitors  of  the  enzyme  could 
serve  as  clinically  useful  antiviral  agents  [2-4].  In  addition,  inhibition  of  UDG 
could  enhance  the  effectiveness  of  current  anticancer  therapies  such  as  5-fluoro- 
uracil  and  methotrexate  that  lead  to  increased  accumulation  of  uracil  in  DNA 
[5,6].  Indeed,  in  yeast  UDG  knockout  strains,  high  levels  of  uracil  in  DNA  leads 
to  cell  cycle  arrest  at  a  G2  checkpoint  [6],  supporting  the  validity  of  this  ap¬ 
proach. 

The  mechanistic  basis  for  the  extraordinary  catalytic  power  of  UDG  has  been 
extensively  investigated  for  the  enzyme  from  Escherichia  coli  and  humans  [1].  A 
key  insight  from  this  work  is  that  the  enzyme  facilitates  a  stepwise  mechanism 
(Fig.  1A),  that  involves  complete  breakage  of  the  A-glycosidic  bond  in  the  first 
transition  state  (TS1),  the  generation  of  a  discrete  oxacarbenium  ion-uracil  anion 
intermediate  [7],  followed  by  attack  of  the  nucleophilic  water  at  Cl'  of  the  inter¬ 
mediate  in  the  second  transition  state  (TS2).  Crystal  structures  have  been  ob¬ 
tained  of  a  reactant  analogue  complex  (U'K,  Fig.  IB)  [8],  a  bimolecular  mimic 
of  TS2  consisting  of  a  cationic  l-aza-2'-deoxyribose  (1-aza-dR)  sugar  and  the  ura¬ 
cil  anion  (I  +  U~,  Fig.  IB),  and  the  reaction  products  of  abasic  DNA  and  uracil 
[9].  The  cationic  1-aza-dR  component  of  the  TS2  mimic  has  been  previously  char¬ 
acterized  as  a  tight  binding  ligand  for  the  UDG-uracil  anion  binary  complex 
( Kd  =  0.5  nM)  [10].  The  high  affinity  of  this  glycosyl  cation  mimic  for  the  binary 
complex  arises  in  large  part  from  favorable  electrostatic  interactions  with  a  con¬ 
served  aspartate  [11],  the  uracil  anion  [11],  and  anionic  DNA  phosphodiester 
groups  [12],  In  contrast,  the  very  similar  (but  neutral)  tetrahydrofuran  abasic  site 
product  mimic  binds  weakly  to  the  enzyme-uracil  anion  complex  (Fig.  IB, 
Kd  >  15  pM)  [10],  establishing  the  importance  of  the  sugar  cation  in  promoting 
tight  binding. 

Although  the  previously  characterized  bipartite  TS2  mimic  is  a  potent  inhibitor  of 
UDG  at  pH  values  in  which  uracil  component  is  anionic  [10,11],  at  neutral  pH  values 
the  affinity  diminishes  greatly,  due  to  the  unfavorable  equilibrium  for  deprotonation 
of  the  base  (pAT™  =  9.8)  [11,13].  To  address  this  shortcoming  we  have  now  synthe¬ 
sized  monopartite  mimics  of  TS1  in  which  the  uracil  base  is  covalently  attached  to 
the  1-azasugar  in  two  ways  (Fig.  IB),  and  we  have  also  explored  several  low  pATa  ura¬ 
cil  analogues  as  improved  coinhibitors  with  1-aza-dR.  We  find  that  one  of  the  new 
TS1  mimics  has  significantly  greater  affinity  than  two  substrate  analogues,  and  that 
the  highest  affinity  TS2  mimic,  consisting  of  1-aza-dR  and  urazole  (pATa  =  5.8),  is  su¬ 
perior  at  neutral  pH  to  the  original  bipartite  construct  using  uracil  (Fig.  IB).  The 
general  merits  and  limitations  of  targeting  TS1  and  TS2  in  stepwise  glycosylase  re¬ 
actions  are  discussed. 
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Fig.  1.  The  reaction  coordinate  of  uracil  DNA  glycosylase  (UDG)  and  reaction  coordinate  mimics.  (A) 
UDG  uses  a  stepwise  mechanism  to  hydrolyze  the  glycosidic  bond  of  deoxyuridine  in  DNA.  The  forma¬ 
tion  and  decay  of  the  oxacarbenium  ion  intermediate  involves  two  high  energy  transition  states  (TS1  and 
TS2)  that  exhibit  unique  geometric  and  electronic  features  that  may  be  mimicked  by  stable  chemical  con¬ 
structs.  (B)  Stable  chemical  mimics  of  UDG  reaction  coordinate  species.  The  2,-fluoro-2'-deoxyuridine 
substrate  analogue  and  the  tetrahydrofuran  abasic  product  analogues  have  been  previously  studied 
[11,15].  The  other  reactant  and  TS  analogues  are  investigated  in  this  work. 


2.  Materials  and  methods 


2.1.  Materials 

As  previously  described,  the  recombinant  UDG  from  E.  coli  strain  B  was  purified 
to  >99%  homogeneity  using  a  T7  polymerase-based  over  expression  system  [14].  The 
concentration  of  the  enzyme  was  determined  using  an  extinction  coefficient  of 
38.5  mM~l  cm"1. 


4 


Y.L.  Jiang  et  al  /  Bioorganic  Chemistry  xxx  (2004)  xxx-xxx 


2.2.  Phosphoramidites 

All  nucleoside  phosphoramidites  were  purchased  from  Applied  Biosystems  or 
Glen  Research  (Sterling,  VA),  except  for  the  2'-P-fluoro-2'-deoxyuridine  phospho- 
ramidite  which  was  synthesized  as  previously  described  [15],  and  the  1-aza-l,  2-dide- 
OXy.4a-carba-D-ribitol  5'-trityl-3'-phosphoramidite,  which  was  synthesized  as 
described  below. 

2.3.  1-Aza-l ,2-dideoxy-4a-carba-n-ribitol  5'-trityl  nucleoside  (2) 

To  a  solution  of  1  (0.6  g,  1.67 mmol)  in  CH2CI2  (10  ml)  [16],  was  added  NEt3 
(0.47  ml,  3.4  mmol)  and  Fmoc-Cl  (0.52 g,  2.0 mmol).  The  reaction  mixture  was 
stirred  under  nitrogen  for  2h,  and  then  purified  by  chromatography  on  silica  with 
ethyl  acetate-hexanes  (1:1,  v/v)  to  give  product  2  (0.78  g)  in  80%  yeild.  'H  NMR 
(CDCI3,  ppm)  5  7.77  (m,  2H);  7.58  (m,  2H);  7.43  (m,  19H);  4.37  (m,  3H);  4.22  (m, 
1H);  3.62  (m,  2H);  3.25  (m,  2H);  3.10  (m,  2H);  and  2.42  (m,  1H). 

2.4.  l-Aza-l,2-dideoxy-4a.-carba-D-ribitol  5' -trityl-3' -phosphoramidite  (3) 

To  a  solution  of  2  (0.226  g,  0.39  mmol)  in  CH2CI2  (10  ml),  was  added  diisopropyl- 
ethylamine  (0.2  ml)  2-cyanoethyl  diisopropylchlorophosphoramidite  (130  pi, 
0.57  mmol).  The  reaction  mixture  was  stirred  under  nitrogen  for  0.5  h,  which  was  pu¬ 
rified  by  chromatography  on  silica  with  ethyl  acetate-hexanes-NEt3  (1:1:0.01,  v/v/v) 
to  give  product  3  (0.22 g)  in  71%  yield.  ’H  NMR  (CDCI3,  ppm)  <5  7.78  (m,  2H);  7.60 
(m,  2H);  7.42  (m,  19H);  4.37  (m,  3H);  4.24  (m,  1H);  3.78  (m,  3H);  3.51  (m,  3H);  3.29 
(m,  1H);  3.10  (m,  3H);  2.60  (m,  3H);  1.14  (m,  12H).  3,P  NMR  (CDC13,  ppm)  5  150.2 
(s).  ESI  calc  for  C48H52N3Na05P  (M  +  Na)  804,  found  804.  This  amidite  is  fairy  un¬ 
stable  in  trace  amounts  of  NEt3,  which  likely  catalyzes  the  cleavage  of  Fmoc  group, 
resulting  in  polymerization  of  the  amidite.  Therefore,  the  amidite  should  be  freshly 
made  for  DNA  synthesis. 

2.5.  Oligonucleotide  synthesis 

The  4  mer  oligonucleotides,  UF,  U'F,  I,  and  <)>  (see  Fig.  IB),  were  synthesized  using 
standard  phosphoramidite  chemistry  with  an  Applied  Biosystems  390  synthesizer.  In 
these  sequences,  UF,  2'-|3-fluoro-2'-deoxyuridine  nucleotide;  U'F,  2'-deoxypseudouri- 
dine  nucleotide;  I,  l-aza-2-dideoxy-4a-carba-D-deoxyribonucleotide;  and  <|>,  tetrahy- 
drofuran  abasic  site  analogue.  During  synthesis  of  oligonucleotide  containing  I,  the 
coupling  time  was  increased  to  10  min.  In  addition,  the  time  for  the  trityl  cleavage 
step  was  increased  to  180s  instead  of  standard  90s.  These  modifications  were  found 
to  increase  the  incorporation  efficiency  at  this  step  from  30  to  80  %.  After  synthesis 
and  deprotection,  the  oligonucleotides  were  purified  by  anion  exchange  HPLC  and 
desalted  by  C-18  reversed  phase  HPLC  (Phenomenex  Aqua  column).  The  correct 
size,  purity,  and  nucleotide  compositions  of  the  final  products  were  assessed  by  an¬ 
alytical  reversed  phase  HPLC  [Phenomenex  Aqua  column  (250  mm  x  10  mm)], 
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MALDI  mass  spectrometry,  and  by  enzymatic  digestion  to  the  constituent  nucleo¬ 
sides  (see  below).  The  concentrations  of  the  4  mer  oligonucleotides  were  determined 
by  UV  absorption  measurements  at  260  nm,  using  the  pair  wise  extinction  coeffi¬ 
cients  for  the  constituent  nucleotides  [17]. 

2.6.  Synthesis  ofTSl  and  TS2  mimics 

As  outlined  in  Fig.  2B,  the  TS1  analogues  with  the  sequences  AU‘AA  and 
AU5AA  were  synthesized  by  derivatizing  the  1-nitrogen  of  AIAA  [12],  using  the  ap¬ 
propriate  uracil  derivative  [18]. 

AU6AA  was  prepared  as  follows:  to  dry  AIAA  (0.9pmol),  were  added  H2O 
(10  pi),  CH3CN  (30  pi),  diisopropylethylamine  (10  pi),  and  6-chloromethyluracil 
(7  pi  of  a  29  mg/ml  solution  in  90:10  MeOH/diisopropylethylamine).  The  clear  solu¬ 
tion  was  stirred  for  3  days  at  room  temperature  and  dried  in  vacuo.  The  final  product 
was  obtained  in  50%  yield  from  the  residue  by  HPLC  using  a  C-18  reverse  phase  col¬ 
umn.  Mass  (ESI)  calc  for  MW  1181,  found  1181. 

AU5AA  was  prepared  as  follows:  to  AIAA  (0.2  pmol  in  20  pi  H2O),  were  added 
uracil  (7  pi  of  a  1.1  mg/ml  solution  H20)  and  HCHO  (9  pi  of  a  0.28  mg/ml  solution 
in  48:48:4  EtOH/CH3CN/NEt3).  The  clear  solution  was  then  dried  in  vacuo.  To  this 
mixture,  CH3OH  (5  pi)  and  CH3CN  (5  pi)  were  added.  The  resulting  mixture  was 


-I  Fmnrn  TrO  -|  _  .Fmoc  c|p.O(CH2)2CN  TrO-i _ Fmoc 

- Fmoci! — - -H'  c  N(iPr)2 _ ^  Ig^N 

0H  NEt3,  CH2CI2, 23  °C,  2  h  0H  DIPEA,  CH2CI2, 23  °C,  30  min  Os  0(CH2)2CN 

,  2  3  N(iPr>* 


Fig.  2.  Synthesis  of  1-aza-deoxyribose  5'-trityl-3'-phosphoramidite  (3)  and  the  synthesis  ofTSl  mimics. 
(A)  3  was  obtained  through  the  sequential  reaction  of  1  with  Fmoc-Cl  and  cyanoethyl-diisopropyl-chlo- 
rophosphoramidite.  (B)  3  was  incorporated  into  the  DNA  sequence  A3AA  (abbreviated  as  I  throughout 
the  text)  using  standard  solid  phase  chemistry.  I  was  used  as  a  common  synthon  for  the  synthesis  of  the 
two  TS1  analogues  U5  and  U6. 
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sealed  and  incubated  at  80-85  °C  for  24  h  and  dried  in  vacuo.  The  final  product  was 
obtained  in  12%  yield  from  the  residue  by  HPLC  using  a  C-18  reverse  phase  column. 
Mass  (ESI)  calc  for  MW  1181,  found  1181. 

2. 7.  Nucleotide  composition  analysis 

The  nucleotide  compositions  of  AU5AA  and  AU6AA  were  confirmed  by  digestion 
with  PI  nuclease  and  alkaline  phosphatase  (both  obtained  from  Roche  Diagnostics) 
followed  by  separation  of  the  constituent  nucleosides  using  reversed  phase  HPLC 
(Phenomenex  C-18  Aqua  column,  5  mm  x  250  mm)  with  monitoring  at  260  nm  and 
isocratic  elution  (7%  CH3CN,  0.1  M  TEAA,  pH  7.0).  The  identity  of  the  peaks 
was  confirmed  by  comparison  with  the  retention  times  of  authentic  nucleoside  stan¬ 
dards  and  the  ratio  of  the  peak  areas  were  consistent  with  the  expected  stoichiome¬ 
tries  and  extinction  coefficients  of  these  tetramers.  The  standards  for  the  21- 
deoxynucleoside  forms  of  U5  and  U6  were  synthesized  as  previously  described  [18]. 
The  structures  of  the  deoxynucleoside  standards  were  ascertained  by  *H-NMR  spec¬ 
troscopy  and  ESI-MS  analyses. 


2.8.  Competitive  inhibition  measurements 

For  measuring  the  binding  of  the  single  stranded  4  mer  UF  substrate  and  <|>  prod¬ 
uct  analogue  DNA  to  the  free  enzyme,  competitive  kinetic  inhibition  measurements 
were  performed  using  the  substrate  ApUpAp  [19].  Conditions  were  chosen  whereby 
[UDG],ot « [inhibitor]  or  [ApUpAp],  and  [ApUpAp]  <Km.  Accordingly,  K,  could  be 
obtained  directly  from  a  plot  of  k/kQ  against  [inhibitor]  as  shown  in  Eq.  (1),  where  k 
is  the  observed  rate  constant  (u/[UDG]tot)  at  a  given  [inhibitor],  and  k0  is  the  ob¬ 
served  rate  constant  in  the  absence  of  the  inhibitor: 

k/ko  =  1/(1  +  [X]/*0.  (1) 

For  these  measurements,  a  sensitive  HPLC  kinetic  assay  for  monitoring  the  for¬ 
mation  of  the  abasic  product  was  employed  [20].  All  experiments  were  performed  at 
25  °C  using  TMN  buffer  at  pH  8  (10  mM  Tris-HCl,  2.5  mM  MgCl2,  and  25  mM 
NaCl).  For  determination  of  the  dissociation  constants  of  the  uracil  analogues  in 
Table  1,  an  alternative  fluorescence-based  competitive  inhibition  kinetic  assay  was 
used  [10], 


2.9.  Binding  of  1-aza-dR  to  the  UDG  binary  complex 

The  dissociation  constants  for  binding  of  I  to  the  UDG-uracil  or  UDG-urazole 
complex  were  determined  by  competition  binding  measurements  in  which  a  2-amino- 
purine  (2AP)  labeled  abasic  analogue  DNA  (4*19)  was  displaced  from  the  EU  or 
EUz  binary  complex  as  previously  described  [10].  The  sequence  of  4>19  has  been  pre¬ 
viously  reported  [21],  and  this  DNA  construct  shows  a  strong  fluorescence  decrease 
when  it  binds  to  the  EU  binary  complex  that  can  be  used  as  a  spectroscopic  signal  in 
competition  binding  measurements.  Measurements  were  performed  at  25  °C  in  TMN 
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buffer  at  pH  7  and  8.  Before  performing  the  competition  binding  experiments,  the  KD 
values  of  4>19  for  the  EU  and  EUz  complexes  were  determined  at  pH  8  and  7  using 
direct  binding  measurements.  In  these  measurements,  the  decrease  in  2AP  fluores¬ 
cence  was  followed  upon  titrating  a  solution  containing  200  nM  <J>19  and  1  mM  ura¬ 
cil  (or  urazole)  with  increasing  concentrations  of  UDG.  This  concentration  of  uracil 
(or  urazole)  is  over  10-fold  greater  than  the  apparent  KD  of  uracil  or  urazole  for  the 
enzyme.  This  insures  that  UDG  is  saturated  with  either  base,  and  that  the  measure¬ 
ments  reflect  binding  of  <(>19  to  the  enzyme-base  binary  complex.  Excitation  was  at 
320  nm  and  emission  spectra  from  340  to  450  nm  were  collected  using  a  Spex  Fluoro- 
max  3  spectrofluorimeter.  The  fluorescence  intensity  (F)  at  370  nm  was  plotted 
against  [EU]tot,  or  [EUz]tot  to  obtain  the  KD  from  Eq.  (2),  where  [X]tot  represents 
either  [EU]tot  or  [EUz],ot. 

F  =  F0-{(F0-  Ff)[4,19]tot/2}{Z>  -  (b2  -  4[X]tot[<|,19]tot)1/2},  (2) 

b  =  KD  +  [X]tot  +  [<|>19]tot. 

To  determine  the  affinity  of  I  for  the  EU  and  EUz  binary  complexes,  titrations 
included  a  saturating  concentration  of  uracil  (1  mM)  so  that  at  the  beginning  of 
the  titration  UDG  was  completely  bound  as  E  •  U  or  E  •  U  •  <j>19.  The  concentrations 
of  <j)19  and  UDG  in  the  individual  experiments  are  reported  in  the  legends  to  Figs.  4 
and  8.  The  dissociation  constants  of  I  for  the  EU  or  EUz  complex  (A^EU)  were  then 
determined  using  the  computer  program  Dynafit  [22]  and  the  equilibria  shown  in 
Eqs.  (3)  and  (4),  employing  the  known  dissociation  constants  of  <|>19  for  the  EU 
and  EUz  complexes  as  determined  from  Eq.  (2)  above: 

jf+.EU 

E  •  U  •  4*19  E  •  U  +  4*19  (3) 

^I.EU 

E  •  U  •  I  ^  E  •  U  + 1  (4) 

2.10.  'H  NMR  spectroscopy 

Samples  (0.5  ml  in  90%  H20  and  10%  D20  for  frequency  lock)  contained  0.3  mM 
UDG,  2mM  uracil,  5-azauracil  or  urazole,  lOmM  NaH2PC>4  (pH  7.5  or  9),  and 
150mM  NaCl.  The  samples  were  placed  in  5  mm  NMR  tubes  (Wilmad  535-PP,  Bue¬ 
na,  NJ)  and  sealed  with  parafilm.  The  experiments  were  collected  at  25  °C  on  a  Var- 
ian  INOVA  500  MHz  spectrometer  using  a  binomial  1-5-10-5-1  pulse  sequence  that 
minimizes  excitation  of  water  [23].  Acquisition  and  processing  parameters  were:  2k 
complex  points,  68  ms  acquisition  time,  and  15  Hz  line  broadening. 

2.11.  Computational  modeling 

The  structural  models  and  electrostatic  potential  surfaces  shown  in  Fig.  9  were 
obtained  with  the  program  Spartan  Pro  (Wavefunction)  using  semi-empirical  meth¬ 
ods  (HF/AM1).  The  truncated  models  included  all  atoms  of  the  base  and  sugar,  but 
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the  5'  and  3'  phosphodiester  groups  were  omitted  and  replaced  with  hydrogen  atoms. 
The  reactant  structure  was  obtained  directly  from  the  coordinates  of  pdb  deposition 
1EMH  without  any  further  optimization  before  calculating  the  electrostatic  potential 
surface.  The  TS1  model  was  obtained  by  constraining  the  sugar  in  a  3'-exo  confor¬ 
mation  calculated  from  KIE  measurements  [7],  and  orienting  the  base  in  the  position 
observed  in  the  structure  of  uracil  and  1-aza-dR  (1QF3),  except  that  the  N1  nitrogen 
was  moved  to  a  distance  of  2.75  A  from  the  anomeric  carbon,  which  corresponds  to  a 
dissociative  transition  state  with  less  than  0.01  bond  order  to  the  leaving  group  [24]. 
The  model  for  TS2  was  obtained  directly  from  the  structure  of  uracil  and  1-aza-dR 
(1QF3)  by  substituting  a  carbon  atom  for  the  1-NHJ  of  the  sugar.  The  geometry  op¬ 
timized  structures  and  electrostatic  potential  surfaces  for  U5  and  U‘  in  Fig.  9B  were 
calculated  using  semi-empirical  methods  (HF/AM1),  and  were  manually  superim¬ 
posed  with  the  structural  models  for  TS1.  The  model  for  bound  urazole  and  1- 
aza-dR  in  Fig.  9C  was  obtained  from  1QF3  by  substituting  the  appropriate  atoms 
of  the  uracil  base,  followed  by  constrained  energy  minimization. 


3.  Results 

3.1.  Synthesis  of  novel  TS1  analogues 

The  characteristics  of  TS1  are  an  elongated  glycosidic  bond,  a  significant  positive 
charge  development  on  the  sugar,  and  negative  charge  development  on  the  uracil 
leaving  group  (Fig.  1A).  To  imitate  these  attributes,  we  synthesized  the  two  TS1 
mimics  shown  in  Fig.  IB.  The  first  mimic  (U6),  has  the  1-nitrogen  linked  to  the  uracil 
base  at  the  6-position  through  a  methylene  bridge.  This  analogue  was  anticipated  to 
imitate  the  elongated  glycosidic  bond  and  the  developing  positive  charge  on  the  su¬ 
gar  in  a  dissociative  transition  state.  In  addition,  since  the  glycosidic  nitrogen  of  the 
uracil  base  of  U6  is  not  involved  in  a  covalent  bond  to  the  sugar,  it  is  free  to  lose  a 
proton  and  generate  an  anion  in  the  active  site,  as  previously  observed  for  uracil 
[13,25],  This  anion  would  be  anticipated  to  provide  stabilization  to  the  glycosyl  cat¬ 
ion  [1 1].  The  second  TS1  mimic  (U5),  is  similar  to  U6  in  that  a  methylene  bridge  con¬ 
nects  the  sugar  and  base  (Fig.  IB).  However,  the  bridge  connects  to  the  5-position  of 
the  uracil,  and  therefore,  the  orientation  of  the  base  heteroatoms  differs  from  U#.  In 
addition,  the  N1  -nitrogen  of  U5  is  not  positioned  as  close  to  the  sugar  as  U6,  and 
therefore  this  position  may  not  ionize  and  provide  stabilization  to  the  glycosyl  cat¬ 
ion. 

The  TS1  analogues  shown  in  Fig.  IB  were  all  synthesized  from  a  common  4  mer 
oligonucleotide  precursor  (I)  that  contained  the  1-azadeoxyribose  moiety  (Figs.  1A 
and  B).  Although  UDG  does  bind  more  tightly  to  longer  oligonucleotides  [10,20], 
the  4  mer  has  been  shown  to  possess  all  of  the  binding  determinants  required  for 
tight  binding  of  the  1-aza-dR  group  [10,12,19].  The  4  mer  was  prepared  using  stan¬ 
dard  solid  phase  phosphoramidite  DNA  chemistry  using  the  commercially  available 
adenosine  nucleotide  3'-phosphoramidite,  and  the  custom  made  l-aza-l,2-dideoxy- 
4a-carba-D-ribitol  5'-trityl-3'-phosphoramidite  (3,  Fig.  1A)  [10,16]. 
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Once  the  I  synthon  was  in  hand  (Fig.  IB),  it  was  fairly  straight  forward  to  deriv- 
atize  its  1 -nitrogen  with  6-chloromethyluracil  or  uracil  in  the  presence  of  formalde¬ 
hyde  to  yield  the  U6  and  U5  TS1  analogues  shown  in  Fig.  2.  These  syntheses  are 
similar  to  those  used  previously  to  obtain  analogous  1-aza-uridine  nucleosides 
[26],  but  this  report  is  the  first  to  generate  these  nucleotides  in  a  DNA  scaffold.  Al¬ 
though  the  yields  are  only  in  the  range  12-30%,  and  the  crude  reaction  products  re¬ 
quired  purification  by  high-performance  liquid  chromatography,  sufficient  material 
was  obtained  to  perform  a  large  number  of  biochemical  experiments. 


3.2.  Relative  binding  affinities  of  substrate,  TS1,  and  bipartite  TS2  analogues 

Two  substrate  analogues,  U'p  and  UF,  were  constructed  for  binding  affinity  com¬ 
parisons  with  the  TS1  and  TS2  analogues  described  below  (Fig.  IB).  Dissociation 
constants  for  the  substrate  and  TS1  mimics  were  determined  by  a  competitive 
inhibition  kinetic  assay  as  previously  described  (Fig.  3)  [12].  U'1’  and  UF  were  found 
to  bind  with  similar  affinities  of  5.5  and  4.8  pM,  respectively,  while  the  TS1  mimic  U5 
bound  10  to  12-fold  more  tightly  (Kn  =  0.5  ±  0.04  pM).  The  other  TS1  mimic,  U®, 
was  found  to  bind  with  similar  affinity  as  the  two  substrate  mimics 
(Kd  =  4.8  ±  0.5  pM),  indicating  that  linking  uracil  via  a  methylene  bridge  to  the 
6-position  is  less  effective  than  to  the  5-carbon  (i.e.,  U5).  These  results  suggest  that 
U5  captures  some  of  the  electronic  and  geometric  features  of  TS1. 

An  assumption  in  ascribing  the  enhanced  binding  affinity  of  U5  to  mimicry  of  TS1 
is  that  the  1 -nitrogen  of  the  sugar  is  protonated.  Previous  NMR  and  pH  studies  have 
established  that  the  p Ka  for  free  1-aza-dR  is  9.5,  and  that  its  nitrogen  is  protonated 
when  bound  to  the  EU  complex  in  the  pH  range  7-9  [11].  Because  direct  measure¬ 
ment  of  the  1 -nitrogen  pA"a  values  in  the  context  of  U5  by  NMR  spectroscopy  is  not 
trivial,  and  we  examined  binding  of  U5  at  pH  6.5  to  assess  whether  the  measurements 
at  pH  8  involved  the  neutral  1-azasugar.  In  contrast  with  the  anticipated  increase  in 
binding  affinity  if  protonation  of  the  sugar  occurred  as  the  pH  was  lowered  from  8.0 


Fig.  3.  Inhibition  of  UDG  by  the  substrate  analogues  U’1’,  UF,  and  the  TS1  analogue  U5.  A  competitive 
kinetic  inhibition  assay  was  used.  The  curves  are  nonlinear  best  fits  to  Eq.  (1). 
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to  6.5,  U5  bound  more  weakly  at  pH  6.5,  with  KD  =  8.3  pM.  The  weaker  binding  at 
pH  6.5  may  be  partially  attributed  to  protonation  of  the  catalytic  aspartate  (p Ka  6.5) 
that  is  positioned  below  the  a  face  of  the  sugar  (Fig.  1A)  [9,27]. 

The  salient  features  of  TS2  are  the  lack  of  covalent  bonding  between  the  sugar 
and  the  anionic  uracil,  and  significant  positive  charge  at  Cl'  (Fig.  1).  Thus  a  good 
mimic  of  TS2  is  the  bipartite  combination  of  uracil  and  I  (Fig.  IB).  Binding  of  I 
to  EU“  complex  can  be  followed  by  competitively  displacing  a  fluorescent  abasic 
DNA  analogue  (<|>,  Fig.  4)  [10],  and  is  pH  dependent  because  interaction  of  this 
glycosyl  cation  mimic  is  enhanced  by  the  negative  charge  on  the  bound  uracil  base 
(pAaEU  =  7.5)  [25].  As  previously  shown,  I  binds  very  tightly  at  pH  8.0  with  a 
£d  =  (5  ±0.4)  xlO”4  pM  (solid  curve,  Fig.  4)  [10].  However  at  pH  7,  where  the 
bound  uracil  is  75%  neutral,  I  binds  300-fold  more  weakly  (A±>  =  0.14  ±  0.03  pM). 
A  bar  chart  is  shown  in  Fig.  5  that  compares  the  relative  affinities  of  the  substrate 
and  TS1  mimics  at  pH  8.0,  and  the  TS2  mimic  at  pH  8.0  and  7.0.  Although  this  com¬ 
parison  clearly  reveals  that  the  bipartite  TS2  mimic  (I  +  U”)  is  superior  at  pH  8.0, 
the  Achilles’  heel  of  this  mimic  is  that  the  pATa  values  for  the  free  and  bound  uracil 
are  too  high  to  generate  the  required  anionic  form  at  neutral  pH,  leading  to  weak 
binding  of  both  U  and  I  under  physiological  conditions.  A  further  discussion  of 
the  merits  and  limitations  of  TS2  mimicry  is  presented  later. 

3.3.  Screening  low  pKa  uracil  analogues  for  improved  coinhibitors 

One  strategy  to  improve  the  potency  of  the  TS2  mimic  at  a  physiological  pH  value 
is  to  find  uracil  analogues  that  have  decreased  p Ka  values  as  compared  to  uracil  itself 


E>1>4  +  I  <  E-U-l  +  ® 


[I] 


Fig.  4.  pH  dependence  of  I  binding  to  the  EU  complex.  A  competitive  displacement  assay  was  used  in 
which  a  fluorescent  abasic  site  analogue  (4>)  was  dissociated  from  the  EU  complex  upon  binding  of  I 
[10].  The  solid  curve  shows  the  previously  obtained  fitted  curve  for  binding  of  I  to  the  EU  complex  at 
pH  8  [10].  The  measurements  at  pH  7  used  1  mM  uracil,  330  nM  UDG,  and  200  nM  <>.  The  curve  was  cal¬ 
culated  using  the  program  Dynafit  [22]. 
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Analogue 

Fig.  5.  Relative  KD  values  (pM)  of  substrate  and  TS1  analogues  at  pH  8.0  (left  graph  axis),  and  the  bipar¬ 
tite  TS2  mimic  at  pH  8.0  and  7.0  (right  graph  axis). 

(p^a(Nl)free  =  9.8,  p^a(Nl)bound  =  6.4)  [13].  With  this  aim,  we  screened  six  uracil  an¬ 
alogues  with  p^a  values  ranging  from  5.8  to  8  for  their  ability  to  competitively  in¬ 
hibit  the  reaction  of  UDG  at  pH  8  and  7  (Table  1).  As  a  point  of  reference,  uracil 
binds  8.2-fold  more  weakly  as  the  pH  is  lowered  in  this  range,  due  to  protonation 
of  the  bound  uracil  [25].  Most  of  the  uracil  analogues  in  Table  1  with  pi£a(Nl)free 
values  of  6.7-8.0  show  1.5-  to  5.9-fold  weaker  binding  as  the  pH  was  lowered  from 
8  to  7,  which  arises  from  protonation  of  both  the  bound  and  free  uracil.  The  mag¬ 
nitude  of  the  binding  decrement  parallels  the  proton  affinities  of  these  analogues  (Ta¬ 
ble  1),  but  it  is  impossible  to  be  more  quantitative  about  the  trends  because  the  p Kz 
values  for  the  bound  uracil  analogues  are  not  known.  One  analogue,  urazole 
(p Ka  =  5.8),  showed  pH  independent  binding  over  the  pH  range  investigated 
(Fig.  6),  and  furthermore,  bound  14-fold  more  tightly  than  uracil  at  pH  7 
(Kd  =  98  pM)  (Fig.  6B,  Table  1). 


Table  1 

p Kd  values  and  dissociation  constants  for  uracil  analogues3 


P*a 

Kd  (mM) 

pH  8 

pH  7 

K*/Kn 

u 

9.8 

0.17  ±0.02 

1.4  ±0.1 

0.12 

5-FU 

8.0 

6.2  ±1.2 

9.6  ±1.6 

0.65 

6-AU 

7.6 

0.39  ±0.01 

1.4  ±0.1 

0.27 

6-CF3-U 

7.4 

0.11  ±0.01 

0.30  ±0.03 

0.37 

CA 

7.2 

0.098  ±0.016 

0.58  ±0.11 

0.17 

5-AU 

6.7 

0.43  ±0.04 

0.98  ±0.02 

0.44 

Uz 

5.8 

0.092  ±0.008 

0.098  ±0.004 

0.94 

aThe  pATa  values  for  U,  5-FU,  6-AU,  6-CF3-U,  CA,  5-AU,  and  Uz  were  obtained  from  [46-52],  re¬ 
spectively. 

bThe  ratio  of  the  dissociation  constants  at  pH  8  and  7. 
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Because  urazole  showed  favorable  binding  properties  to  the  free  enzyme,  we  in¬ 
vestigated  its  binding  interactions  in  more  detail  using  NMR  spectroscopy.  Previous 
NMR  studies  of  uracil  binding  to  UDG  have  revealed  the  presence  of  a  very  down- 
field  shifted  proton  resonance  that  arises  from  a  hydrogen  bond  between  uracil  02 
and  the  imidazole  NH  of  Hisl87  in  the  active  site  (Fig.  1A)  [11,13,28].  This  hydrogen 
bond  stabilizes  the  uracil  anion  in  the  active  site  by  5  kcal/mol  and  plays  an  indirect 
but  essential  role  in  binding  of  I  by  facilitating  formation  of  the  uracil  anion  compo¬ 
nent  of  the  electrostatic  sandwich  that  cradles  the  sugar  cation  [11],  To  directly  test 
whether  the  urazole  anion  persists  at  neutral  pH  in  the  UDG  active  site,  we  collected 
NMR  spectra  of  UDG  in  the  presence  of  uracil,  5-azauracil,  and  urazole  at  pH  9.0 
and  7.5  (Figs.  7A  and  B).  At  pH  9.0,  uracil  and  urazole  both  show  the  expected  res¬ 
onance  at  14.6  and  14.8  ppm,  but  5-azauracil  does  not.  The  absence  of  the  proton 
resonance  for  5-azauracil  is  consistent  with  its  weaker  binding  to  the  enzyme,  and 
suggests  that  this  may  be  due  to  its  improper  positioning  in  the  active  site.  When 
the  pH  is  lowered  to  7.5,  only  the  downfield  resonance  of  urazole  remains,  while  that 
for  uracil  is  broadened  beyond  the  limits  of  detection.  This  NMR  result  strongly  sup¬ 
ports  the  proposal  that  the  pH  independent  binding  of  urazole,  and  its  enhanced  af- 


Fig.  6.  Competitive  inhibition  of  the  reaction  of  UDG  with  the  substrate  AUAp  (1  pM)  by  uracil  and 
urazole  at  (A)  pH  8  and  (B)  pH  7. 
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A  pH  9.0 


U 

AzU 


15.8  15.4  15.0  14.6  14.2  13.8  13.4  ppm 


B  pH  7.5 


u 


AzU 


15.8  18.4  15.0  14.6  14.2  13.8  13.4  ppm 


Fig.  7.  *H  NMR  spectra  of  uracil  (U),  5-azauracil  (AzU),  and  urazole  (Uz)  bound  to  UDG  at  (A)  pH  9.0 
and  (B)  pH  7.5.  The  downfield  resonance  arises  from  the  interaction  of  the  NeH  of  His  187  with  the  02 
anion  of  U  and  Uz  (see  text). 


EHJz>j)  +  I  <  -  E’Uz'l  +  4 


Fig.  8.  Binding  of  I  to  the  EUz  complex  at  pH  8  and  7.  The  measurements  at  pH  8  used  1.4pM  UDG, 
200  nM  <|>,  and  1  mM  urazole  (Uz).  The  measurements  at  pH  7  used  330  nM  UDG,  200  nM  4>,  and  1  mM 
urazole  (Uz).  The  curves  were  calculated  with  the  program  Dynafit  [22]. 

finity  as  compared  to  uracil,  arises  from  its  reduced  p Ka,  allowing  it  to  be  negatively 
charged  at  neutral  pH. 

Because  of  the  enhanced  binding  of  urazole  at  neutral  pH,  we  investigated  this 
analogue  as  a  coinhibitor  with  I  (Fig.  8).  Using  the  competitive  displacement 
fluorescence  assay,  I  was  found  to  bind  to  the  EUz-  complex  with  an  affinity  of 
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14  ±  4  and  49  ±  5  nM  at  pH  8  and  7,  respectively.  Thus  at  a  physiological  pH  7.4, 1 
binds  to  the  EUz~  complex  with  an  affinity  of  about  30  nM  as  estimated  from  inter¬ 
polation  between  these  measured  values.  This  affinity  is  170-fold  greater  than  the 
substrate  analogues  bind  to  the  free  enzyme. 


4.  Discussion 

4.1.  Targeting  TS1  in  stepwise  glycosylase  reactions 

Glycosylase  reactions  have  long  been  amenable  to  the  design  of  inhibitors  that  mi¬ 
mic  features  of  high  energy  structures  along  the  reaction  coordinate,  thereby  captur¬ 
ing  a  portion  of  the  enzyme’s  strong  binding  energy  for  these  species  [29-35].  In  the 
case  of  a  stepwise  glycosylase  reaction,  where  a  discrete  oxacarbenium  ion  interme¬ 
diate  is  formed,  it  is  possible  to  envision  the  design  of  molecules  that  mimic  TS1,  the 
intermediate,  or  TS2.  In  general,  which  of  these  mimics  becomes  the  most  useful  in¬ 
hibitor  will  depend  on  how  closely  the  unique  charge  and  geometric  properties  of 
these  reaction  coordinate  species  are  reproduced  by  the  inhibitor.  In  addition,  to  be¬ 
come  useful  inhibitors,  bipartite  intermediate  and  TS2  mimics  must  also  possess  suf¬ 
ficient  binding  energy  to  overcome  the  unfavorable  entropic  penalties  associated 
with  binding  of  the  two  molecular  components  from  solution.  A  transition  state 
mimicry  strategy  employing  both  1'  and  4'-azasugars  has  been  used  previously  to  po¬ 
tently  inhibit  several  DNA  glycosylases  [36,37]. 

Mimics  of  TS1  in  which  the  sugar  and  base  are  covalently  tethered  must  reproduce 
the  elongated  bond  lengths  of  this  dissociative  transition  state  [24,38],  as  well  as  the 
charge  distributions  on  the  sugar  and  base.  In  addition,  the  linker  must  allow  favorable 
positioning  of  the  leaving  group  and  sugar  relative  to  enzyme  groups  that  form  stabi¬ 
lizing  interactions  with  these  parts  of  the  substrate.  Using  the  classic  pseudo-thermo¬ 
dynamic  framework  developed  by  Radzicka  and  Wolfenden  [39],  the  maximum 
binding  affinity  for  a  perfect  TS1  mimic  would  be  equal  to  KD  =  knon/(kcat /Km),  where 
knon  is  the  rate  of  the  uncatalyzed  glycosidic  bond  cleavage  reaction.  In  the  case  of  UDG 
and  its  substrate  AUAA,A:non/(fcCat/^m)  =  10-10  s~'/(3  x  106M-1  s"')~  10_17M[1,10]. 
Accordingly,  a  perfect  TS1  mimic  would  be  expected  to  bind  almost  twelve-orders  of 
magnitude  more  tightly  than  the  Uy  and  UF  substrate  analogues  (Fig.  IB).  Although 
the  TS1  mimics  studied  here  bind  as  much  as  12-fold  more  tightly  than  the  substrate 
analogues,  they  capture  only  a  tiny  fraction  of  the  theoretical  interaction  energy  ex¬ 
pected  from  a  perfect  TS1  mimic. 

Previous  structural  and  mechanistic  studies  provide  an  informative  basis  for  un¬ 
derstanding  the  comparatively  weak  binding  affinity  of  the  TS1  mimics.  In  the  crystal 
structure  of  UDG  bound  to  DNA  containing  the  substrate  analogue  deoxypseudo- 
uridine  (U'1’),  the  nucleotide  base  and  sugar  are  oriented  in  a  highly  distorted  confor¬ 
mation  in  which  the  base  is  nearly  coplanar  with  the  sugar  ring,  and  the  C-C 
glycosidic  bond  is  elongated  (1.55  A,  Fig.  9A)  [8].  This  ground-state  conformation 
appears  to  be  significantly  strained,  and  likely  represents  a  high  energy  conformation 
approaching  that  of  TS1,  which  is  highly  dissociative  [7].  A  reasonable  structural 
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A  Reactant  TS1  TS2 


b  TS1  mimics 


C  TS2  mimic 


Fig.  9.  Structural  and  electrostatic  potential  models  for  reactant,  TS1,  TS2,  and  chemical  mimics.  (A) 
Structural  models  for  the  bound  substrate  analogue  li4*  and  the  two  transition  states  of  the  stepwise  re¬ 
action  catalyzed  by  UDG.  (B)  Geometry  optimized  structural  models  for  Us  and  V6  superimposed  on 
the  model  for  TS1.  The  electrostatic  potentials  are  plotted  on  the  van  der  Waals  surfaces,  which  are  shown 
to  the  right  of  each  model.  (C)  Geometry  optimized  structural  model  for  the  bipartite  TS2  mimic  of  uraz- 
ole  and  I  superimposed  on  the  model  for  TS2.  The  electrostatic  potential  surface  of  Uz  4- 1  is  shown  to  the 
right.  The  details  of  how  these  models  were  calculated  are  described  in  Section  2. 


model  for  TS1  is  shown  in  Fig.  9 A,  which  was  obtained  from  three  pieces  of  infor¬ 
mation  (see  Section  2):  (i)  the  crystallographic  coordinates  of  UT  [8],  (ii)  the  results 
from  KIE  studies  that  provide  the  sugar  pucker  and  Cl'-Nl  bond  order  in  TS1 
(<0.01)  [7],  and  the  position  of  the  uracil  base  in  the  crystal  structure  of  the  bipartite 
TS2  mimic  of  uracil  and  1-aza-dR  [40].  A  model  of  TS1  obtained  using  this  informa¬ 
tion  is  superimposed  with  energy  minimized  models  of  U6  and  U5  in  Fig.  9B.  As 
compared  to  the  glycosidic  bond  of  the  substrate  analogue  U*,  the  methylene 
bridges  of  U6  and  U5  increase  the  linear  distance  between  the  base  and  sugar  from 
1.55  to  about  2.6  A.  However  as  shown  in  Fig.  9B,  the  constraints  imposed  by  the 
methylene  bridge  groups  do  not  allow  an  orientation  of  the  base  and  sugar  that  pre¬ 
cisely  mimics  that  of  TS1. 

In  addition  to  the  relatively  poor  structural  mimicry  by  U6  and  U5,  these  ana¬ 
logues  also  poorly  match  the  charge  distribution  on  the  sugar  and  base  in  TS1. 
This  can  be  appreciated  from  comparison  of  the  electrostatic  potential  surfaces  of 
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these  TS1  mimics  with  that  of  the  TS1  model  (Figs.  9A  and  B).  For  instance,  U5  and 
U6  both  mimic  the  positive  charge  character  at  Cl',  but  they  poorly  match  the  neg¬ 
ative  potential  on  the  uracil  leaving  group,  because  the  electrons  of  their  intact  gly- 
cosidic  bonds  are  not  able  to  delocalize  onto  the  02  atom  of  the  base.  A  potent 
purine  based  TS1  inhibitor  (4)  analogous  to  Us  and  U6  has  been  synthesized  for 
the  enzyme  purine  nucleoside  phosphorylase  (PNP)  from  bovine  and  Mycobacterium 
tuberculosis  (Mt)  [34,41], 


The  MtPNP  follows  a  dissociative  but  concerted  mechanism  and  the  methylene 
bridge  analogue  (4)  presumably  mimics  the  features  of  this  transition-state 
(Kd  =  24  pM)  [34].  In  contrast,  the  bovine  PNP  shows  an  earlier  more  associative 
transition-state,  and  accordingly,  (4)  bound  much  more  weakly  [34,42].  For  UDG, 
neither  type  of  TS1  inhibitor  is  very  effective  because  of  the  highly  unusual  orienta¬ 
tion  of  the  sugar  and  base,  and  the  optimization  of  the  active  site  towards  stabiliza¬ 
tion  of  the  electronic  features  of  the  oxacarbenium  ion  intermediate  and  TS2. 

4.2.  The  merits  of  targeting  TS2 

The  properties  of  TS2  are  a  fully  dissociated  uracil  base,  which  has  accumulated  a 
full  negative  charge,  and  a  glycosyl  cation  that  will  have  developed  some  bonding  to 
the  incoming  water  nucleophile.  The  recent  crystal  structure  of  I  and  the  uracil  anion 
bound  to  human  UDG  shows  that  I  assumes  a  distorted  conformation  in  which  the 
imino  nitrogen  is  displaced  from  the  plane  of  the  sugar  ring  towards  the  attacking 
water  [40].  An  implication  of  this  structure  is  that  the  planar  conformation  of  the 
oxacarbenium  ion  intermediate  is  broken  in  TS2,  and  that  the  electrophilic  anomeric 
carbon  migrates  to  meet  the  water  nucleophile.  Thus,  although  I  was  originally 
thought  to  mimic  the  planar  oxacarbenium  ion  intermediate,  its  nonplanar  confor¬ 
mation,  and  its  charge  localization  on  the  1'  position,  makes  it  a  much  better  mimic 
of  TS2.  The  structure  and  corresponding  electrostatic  model  of  the  bipartite  TS2  mi¬ 
mic  consisting  of  urazole  and  I  is  shown  in  Fig.  9  for  comparison  with  a  model  of 
TS2  (Fig.  9A),  which  was  constructed  from  the  crystallographic  coordinates  of 
1Q3F.  The  excellent  geometric  and  electrostatic  match  between  the  bipartite  inhibi¬ 
tor  and  TS2  is  much  better  than  that  observed  for  any  of  the  TS1  mimics  (see  above). 

The  development  of  bipartite  TS2  inhibitors  for  glycosylase  reactions  is  a  nontra- 
ditional  approach,  because  a  single  substrate  molecule  is  deconstructed  into  two 
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parts,  each  of  which  mimics  the  presumed  features  of  a  bimolecular  transition  state. 
This  contrasts  with  the  bisubstrate  analogue  strategy,  where  two  reactants  are  cova¬ 
lently  tethered  to  capture  the  binding  energy  of  the  two  halves,  thereby  bypassing  the 
entropic  penalty  for  binding  two  molecules  from  solution  [43,44].  The  inherent  dif¬ 
ficulty  with  this  bimolecular  approach  is  that  the  binding  of  each  component  is  com¬ 
paratively  weak  in  the  absence  of  the  other  component.  Thus  in  general,  the 
individual  binding  energy  of  at  least  one  of  the  parts  must  be  sufficiently  large  to  sig¬ 
nificantly  populate  the  enzyme  under  physiological  conditions.  Under  such  condi¬ 
tions,  binding  of  the  first  component  creates  the  binding  environment  that  allows 
tight  binding  of  the  second  component,  and  creation  of  the  high  affinity  bimolecular 
inhibitor  complex.  In  practical  terms,  this  requires  that  the  intracellular  concentra¬ 
tion  of  at  least  one  of  the  component  parts  be  high  relative  to  its  KD.  Of  course,  a 
monopartite  TS2  mimic  could  provide  a  significant  entropic  advantage  over  the  bi¬ 
partite  mimic  if  the  correct  structural  features  could  be  captured. 

In  principle,  the  above  thermodynamic  problem  for  a  bipartite  inhibitor  may  be 
overcome  if  one  of  the  components  is  also  the  last  reaction  product  to  be  released 
(such  as  the  uracil  anion  in  this  case)  [10].  For  such  an  ordered  product  release  mech¬ 
anism,  the  second  component  (1-aza-dR  in  this  case)  will  bind  avidly  to  the  enzyme- 
product  complex,  and  the  inhibition  mechanism  will  be  noncompetitive  with  respect 
to  substrate.  For  UDG,  in  vitro  studies  at  pH  8.0  have  shown  that  80  nM  I  provided 
95%  inhibition  of  UDG  in  the  absence  of  any  added  uracil  coinhibitor,  and  that  this 
strong  inhibition  arose  from  I  binding  tightly  to  the  EU  product  complex  [10].  Of 
course,  both  U  and  I  also  bind  to  the  free  enzyme,  but  these  inhibition  pathways  will 
not  be  significant  given  the  weak  affinities  of  the  individual  parts  for  the  free  enzyme 
at  physiological  pH.  A  noncompetitive  mode  of  inhibition  is  especially  advantageous 
for  targeting  a  DNA  repair  enzyme  that  not  only  binds  to  damaged  target  sites  in 
DNA,  but  also  binds  nonspecifically  to  nontarget  DNA,  which  is  present  at  high 
concentrations  in  cells. 


5.  Conclusion 

We  have  explored  chemical  analogues  that  mimic  the  features  of  several  species 
along  the  stepwise  enzymatic  reaction  coordinate  of  UDG.  In  general,  these  findings 
show  that  bipartite  TS2  mimics  can  serve  as  effective  inhibitors  for  glycosylase  reac¬ 
tions  that  proceed  by  stepwise  mechanisms.  In  favorable  cases,  the  bipartite  ap¬ 
proach  can  offer  advantages  over  unimolecular  TS1  inhibitors  because  stable 
chemical  analogues  of  TS2  may  be  better  able  to  match  the  geometries  and  charge 
distributions  of  the  cationic  sugar  and  leaving  group.  Furthermore,  we  suggest  that 
this  strategy  is  likely  to  be  most  successful  when  the  following  conditions  are  met:  (i) 
when  one  of  the  inhibitor  components  is  the  anionic  product  leaving  group  of  the 
reaction,  (ii)  when  the  free  and  bound  forms  of  the  leaving  group  product  have 
the  correct  anionic  ionization  state  at  physiological  pH,  and  (iii)  when  the  intracel¬ 
lular  concentration  of  the  leaving  group  product  is  comparable  to  or  greater  than  its 
Kd.  In  addition,  the  inhibition  is  noncompetitive  with  respect  to  substrate,  which  is 
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desirable  when  high  concentrations  of  substrate  are  present  in  the  cellular  environ¬ 
ment.  Several  reactions  that  potentially  meet  these  criteria  are  those  catalyzed  by  su- 
gar-NDP  hydrolases  or  transferases  [32,45],  nucleoside  phosphorylases  [42],  and 
phosphoribosyl  transferases  [42]. 
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abstract:  To  efficiently  maintain  their  genomic  integrity,  DNA  repair  glycosylases  must  exhibit  high 
catalytic  specificity  for  their  cognate  damaged  bases  using  an  extrahelical  recognition  mechanism.  One 
possible  contribution  to  specificity  is  the  weak  base  pairing  and  inherent  instability  of  damaged  sites 
which  may  lead  to  increased  extrahelicity  of  the  damaged  base  and  enhanced  recognition  of  these  sites. 
This  model  predicts  that  the  binding  affinity  of  the  enzyme  should  increase  as  the  thermodynamic  stability 
of  the  lesion  base  pair  decreases,  because  less  work  is  required  to  extrude  the  base  into  its  active  site.  We 
have  tested  this  hypothesis  with  uracil  DNA  glycosylase  (UDG)  by  constructing  a  series  of  DNA  duplexes 
containing  a  single  uracil  (U)  opposite  a  variety  of  bases  (X)  that  formed  from  zero  to  three  hydrogen 
bonds  with  U.  Linear  free  energy  (LFE)  relationships  were  observed  that  correlated  UDG  binding  affinity 
with  the  entropy  and  enthalpy  of  duplex  melting,  and  the  dynamic  accessibility  of  the  damaged  site  to 
chemical  oxidation.  These  LFEs  indicate  that  the  increased  conformational  freedom  of  the  damaged  site 
brought  about  by  enthalpic  destabilization  of  the  base  pair  promotes  the  formation  of  extrahelical  states 
that  enhance  specific  recognition  by  as  much  as  3000-fold.  However,  given  the  small  stability  differences 
between  normal  base  pairs  and  U*A  or  U*G  base  pairs,  relative  base  pair  stability  contributes  little  to  the 
>  106-fold  discrimination  of  UDG  for  uracil  sites  in  cellular  DNA.  In  contrast,  the  intrinsic  instability  of 
other  more  egregious  DNA  lesions  may  contribute  significantly  to  the  specificity  of  other  DNA  repair 
enzymes  that  bind  to  extrahelical  bases. 


The  genetic  information  of  a  cell  can  be  irreversibly  altered 
through  the  chemical  modification  of  nucleotide  bases  (7). 
To  combat  these  mutagenic  effects,  organisms  have  evolved 
a  two-tiered  base  excision  repair  (BER)  pathway  that  handles 
a  wide  array  of  base  lesions  (2).  In  the  first  stage,  a  highly 
specific  DNA  repair  glycosylase  excises  the  damaged  base 
from  the  DNA  ( 3 ,  4),  producing  an  abasic  site.  This 
intermediate  is  then  processed  by  the  sequential  action  of 
several  repair  enzymes  that  ultimately  restore  the  site  to  its 
original  state  (2).  Since  damaged  site  specificity  resides  solely 
with  the  DNA  glycosylase,  these  enzymes  must  possess 
extraordinarily  high  catalytic  specificities  (4,  5).  In  the 
absence  of  such  specificity,  undamaged  bases  would  be 
randomly  excised  from  DNA,  leading  to  undesirable  abasic 
sites  and  genetic  instability  (6-10). 

In  general,  enzymatic  specificity  results  from  the  extra¬ 
ordinary  structure  of  enzyme  active  sites  that  disfavors 
formation  of  catalytically  productive  interactions  with  non¬ 
substrate  molecules,  and  strongly  favors  such  interactions 
with  the  true  substrates.  As  an  essential  part  of  their 
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recognition  mechanisms,  all  DNA  glycosylases  extrude  their 
damaged  bases  from  the  DNA  double  helix  in  a  process 
known  as  base  flipping,  thereby  placing  it  extrahelically  into 
their  active  sites  where  specific  interactions  with  the  damaged 
base  can  be  formed  (77,  72).  As  part  of  the  energetic  cost  of 
base  flipping,  the  hydrogen  bonds  and  stacking  interactions 
of  the  base  pair  must  be  disrupted.  Thus,  a  prediction  is  that 
DNA  glycosylases  should  bind  more  tightly  to  damaged  sites 
with  disrupted  base  pairing  because  it  requires  less  binding 
energy  to  flip  the  damaged  base  from  the  destabilized  site 
(4,  13—22).  Such  a  thermodynamic  mechanism  is  quite 
general,  and  would  apply  even  for  DNA  glycosylases  that 
interact  with  the  base  that  opposes  the  damaged  base  (23- 
25). 

How  much  does  the  intrinsic  thermodynamic  stability  of 
the  damaged  base  pair  contribute  to  specific  damaged  site 
binding  by  DNA  glycosylases?  We  have  investigated  this 
question  using  the  enzyme  uracil  DNA  glycosylase  (UDG),1 
which  removes  uracil  from  U*G  and  U*A  base  pairs  in  duplex 
DNA  (74, 26, 27).  The  approach  was  to  measure  the  binding 
affinity  of  UDG  for  a  series  of  DNA  duplexes,  in  which  the 
number  of  hydrogen  bonds  (n)  between  uracil  and  its 
opposing  base  (X)  were  systematically  varied  (Figure  1A). 


1  Abbreviations:  UDG,  uracil  DNA  glycosylase;  UF,  2'-/?-fluoro-2'- 
deoxyuridine;  <j>,  abasic  site;  D,  2,6-diaminopurine;  M,  4-methylindole; 
N,  nebularine;  DSC,  differential  scanning  calorimetry;  \p,  pseudo¬ 
dihedral  angle. 
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Figure  1:  (A)  Specific  binding  of  UDG  (E)  to  a  DNA  duplex  containing  deoxyuridine  (U)  opposite  a  purine  analogue  (X).  In  this  study, 
the  number  of  hydrogen  bonds  (n)  in  the  U*X„  base  pair  varies  from  zero  to  three.  Since  binding  of  uracil  requires  breaking  of  the  base  pair 
hydrogen  bonds,  the  difference  in  binding  free  energy  (AAGbind)  between  a  DNA  duplex  with  n  hydrogen  bonds  and  a  reference  duplex 
with  three  hydrogen  bonds  (U*X3)  should  in  part  reflect  the  reduced  energetic  cost  of  breaking  the  hydrogen  bonds  in  the  base  pair.  (B) 
Dissociation  of  a  duplex  with  a  U*X„  base  pair  into  two  single  strands.  The  difference  in  the  duplex  melting  free  energy  (AAGmeit)  between 
a  U-X*  duplex  and  the  reference  duplex  (U*X3)  will  also  reflect  the  energetic  effects  of  a  destabilized  base  pair. 


Then,  these  binding  affinities  were  correlated  with  rigorous 
measurements  of  the  thermodynamic  stabilities  of  these 
duplexes  (Figure  IB).  UDG  was  selected  for  this  study 
because,  unlike  many  DNA  glycosylases  (23—25),  it  does 
not  make  any  direct  contacts  with  the  base  opposite  uracil 
(11).  Thus,  the  observed  changes  in  binding  affinity  can  be 
largely  attributed  to  the  relative  stability  of  the  U*X„  base 
pairs  in  the  free  duplex  DNA  and  not  differential  interactions 
of  the  enzyme  with  the  various  opposing  bases.  These 
quantitative  free  energy  correlations  provide  the  first  direct 
evidence  that  promotion  of  extrahelical  conformations  by 
enthalpic  destabilization  of  a  damaged  site  can  indeed 
enhance  the  specific  binding  of  a  DNA  repair  enzyme. 

EXPERIMENTAL  PROCEDURES2 

Materials.  The  2'-deoxynucleoside  phosphoramidites,  CPG 
supports,  and  DNA  synthesis  reagents  were  purchased  from 
Glen  Research  (Sterling,  VA),  except  for  2'-/7-fluoro-2'- 
deoxyuridine  (UF),  which  was  synthesized  as  described 
previously  (5,  28).  The  oligonucleotides  were  synthesized 
using  standard  phosphoramidite  chemistry  on  an  Applied 
Biosystems  392  synthesizer.  The  oligonucleotides  were 
purified  by  anion  exchange  HPLC  (Zorbax),  followed  by 
C-18  reversed  phase  HPLC  (Phenomenex  Aqua  column). 
Fractions  containing  pure  oligonucleotide  were  concentrated 
to  dryness  in  vacuo ,  redissolved  in  MilliQ  water,  and  stored 
at  —20  °C.  The  purity  of  the  oligonucleotides  was  assessed 
by  matrix-assisted  laser  desorption  mass  spectroscopy  and 
denaturing  polyacrylamide  gel  electrophoresis.  The  concen¬ 
tration  of  each  oligonucleotide  was  determined  using  its 
extinction  coefficient  at  260  nm  (29).  DNA  duplexes  were 
hybridized  in  10  mM  Tris-HCl  (pH  8.0)  and  25  mM  NaCl 
as  described  previously  (5).  The  purification  of  Escherichia 
coli  UDG  has  been  described  previously  (30). 

Kd  Measurements.  The  KD  values  for  binding  of  the  UF- 
X„  duplexes  to  UDG  were  measured  essentially  as  described 
using  a  kinetic  competitive  inhibition  HPLC  assay  under 
conditions  where  the  apparent  Kx  is  equal  to  the  Kd  value 
(i.e.,  [S]  «:  Km,  where  S  exhibits  rapid  equilibrium  binding) 
(31).  The  only  modification  was  that  the  abasic  product 
.  (Ad>Ap)  and  reactant  (AUAp)  were  separated  using  isocratic 


2  Certain  commercial  materials,  instruments,  and  equipment  are 
identified  herein  to  specify  the  experimental  procedure  as  completely 
as  possible.  In  no  case  does  such  identification  imply  a  recommendation 
or  endorsement  by  the  National  Institute  of  Standards  and  Technology, 
nor  does  it  imply  that  the  material,  instruments,  or  equipment  identified 
is  necessarily  the  best  available  for  the  purpose. 


conditions  with  9.5%  CH3CN  and  0.1  M  triethyl  ammonium 
acetate.  Reaction  mixtures  (35  / iL )  containing  10  mM  Tris- 
HCl  (pH  8.0),  60  mM  NaCl,  12.5  ^g/mL  BSA,  1  pM  AUAp, 
0.5  nM  UDG,  and  a  variable  amount  of  the  UF*X„  duplex 
were  incubated  at  room  temperature  for  5  min.  The  KD  for 
each  duplex  was  determined  by  fitting  to  eq  1 

k'K  =  1/(1  +  [UF-XJ/KD)  (1) 

where  k\  is  the  inhibited  rate  and  k0  is  the  rate  in  the  absence 
of  competitor  DNA.  For  the  tightest  binding  duplex  (UF- 
M),  eq  1  was  modified  to  take  into  account  inhibitor  depletion 
(32). 

The  differences  in  binding  free  energies  relative  to  the 
duplex  with  three  hydrogen  bonds  (UF,D)  were  calculated 
from  the  measured  KD  values  using  eq  2. 

AAGbind  =  RT  ln[XD(UF-X„)/KD(UF-D)]  (2) 

Fluorescence  Spectroscopy.  To  ascertain  that  all  of  the 
DNA  duplexes  attained  the  same  bound  state,  tryptophan 
fluorescence  measurements  of  free  and  DNA-bound  UDG 
were  performed.  Samples  (497  pV)  containing  10  mM  Tris- 
HCl  (pH  8.0),  60  mM  NaCl,  and  300  nM  UDG  were 
incubated  for  3  min  at  25  °C  in  a  10  mm  quartz  cuvette, 
and  a  fluorescence  emission  spectrum  was  recorded  in  the 
range  of  325-425  nm  on  a  SPEX  FluoroMax-3  fluorimeter 
(Aex  =  295  nm).  UF*X„  duplex  DNA  (2.5  pL)  was  then  added 
to  the  UDG  solution  to  give  a  final  DNA  concentration  of 
500  nM.  The  reaction  mixture  was  magnetically  stirred  and 
incubated  for  3  min  at  25  °C,  before  the  fluorescence 
emission  spectrum  of  the  UDG-DNA  complex  was  recorded. 
The  tryptophan  fluorescence  intensities  of  free  and  DNA- 
bound  UDG  at  333  nm  were  measured,  and  the  raw  values 
were  then  normalized  for  the  fraction  of  UDG  bound  to  each 
DNA  analogue  before  the  ratio  (i^33bound/F^33free)  was  calcu¬ 
lated. 

Differential  Scanning  Calorimetry.  DSC  measurements  of 
duplex  strand  melting  were  taken  using  a  VP-DSC  micro¬ 
calorimeter  from  Microcal,  Inc.  (Northampton,  MA)  es¬ 
sentially  as  described  previously  (33).  The  DNA  solutions 
had  a  concentration  of  20  pM  in  DNA  duplex  with  10  mM 
Na2HP04  (pH  7.5)  and  60  mM  NaCl.  Samples  were 
equilibrated  at  20  °C  for  15  min  and  scanned  up  to  95  °C  at 
a  preset  scan  rate  of  60  °C/h.  The  transition  peak  areas  were 
measured  using  the  EXAM  software  program  (34),  and  the 
transition  peak  areas  were  divided  by  the  DNA  duplex 
concentration  to  provide  the  transition  enthalpies.  Transition 
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Figure  2:  Structures  of  UF-X„  base  pairs:  UF*D  (54,  55),  UF*G 
(5(5),  UF,N  (57,  55),  UF-<F  (5P-<57),  and  UF-M  ((52). 


entropies  were  determined  from  plots  of  CpfT  against  T,  by 
dividing  the  transition  peak  area  by  the  DNA  duplex 
concentration.  Procedures  for  buffer  baseline  corrections  and 
accounting  for  pre-  and  post-transition  baselines  have  been 
described  previously  (33). 

KMnOj  Oxidation  Measurements.  Because  of  the  low 
reactivity  of  UF  to  oxidation  (55),  it  was  replaced  with 
thymine  (T)  in  the  oligonucleotides  used  in  this  study.  To 
reaction  mixtures  (20  fiL)  containing  10  mM  Tris-HCl  (pH 
8.0),  60  mM  NaCI,  and  either  100  nM  single-stranded  5'- 
[32P]T  or  5'-[32P]T*Xw  duplex  was  added  2.5  mM  KMn04. 
After  the  sample  had  been  incubated  for  3  min  at  room 
temperature,  the  reaction  was  halted  by  the  addition  of  20 
fiL  of  a  stop  solution  containing  1.5  M  sodium  acetate,  1  M 
2-mercaptoethanol,  and  200  pig/mL  tRNA.  The  samples  were 
processed,  imaged,  and  quantified  as  described  previously 
(32). 

RESULTS 

Binding  of  UDG  to  Destabilized  Damaged  Sites.  A  series 
of  15mer  duplexes  were  constructed  in  which  2'-/?-fluoro- 
2'-deoxyuridine  (UF),  a  nonhydrolyzable  uracil  analogue  (5), 
was  placed  opposite  a  series  of  bases  (X)  which  form  zero 
to  three  hydrogen  bonds  with  UF  (Figure  2).  The  affinity  of 
UDG  for  each  UF*X  duplex  was  measured  by  a  competitive 
inhibition  kinetic  assay  in  which  a  3mer  substrate  (AUAp) 
is  separated  from  the  abasic  product  (A4>Ap)  using  reverse 
phase  HPLC  (Figure  3A)  (31).  Representative  inhibition  data 
are  shown  in  Figure  3B  for  the  duplex  that  contains  a 
uracil *4-methyIindole  base  pair  (UF,M).  The  Kb  values  for 


Figure  3:  Determination  of  the  binding  affinity  of  specific  DNA 
duplexes  using  a  competitive  inhibition  assay.  (A)  UDG  HPLC 
activity  assay.  The  substrate  (AUAp)  and  product  A<FAp  are 
indicated.  (B)  Inhibition  of  UDG  by  a  U*M  duplex  (KD  =  0.24  ± 
0.03  nM). 
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Figure  4:  (A)  Binding  affinities  of  UDG  for  UF*X„  DNA  duplexes. 
The  Kb  values  (nanomolar)  are  shown  above  the  individual  bars. 
(B)  Quench  of  UDG  tryptophan  fluorescence  upon  binding  of  the 
UF*X„  duplex.  The  degree  of  UDG  quenching  was  normalized  to 
reflect  the  quenching  when  UDG  is  saturated  with  DNA  (see 
Experimental  Procedures). 

the  six  duplexes  showed  an  incremental  decrease  as  the 
number  of  hydrogen  bonds  was  ablated  (Figure  4A).  The 
tightest  observed  binding  affinity  was  for  the  UF-M  construct 
that  has  no  hydrogen  bonds  (Kb  =  0.25  nM).  This  affinity 
is  3000-fold  tighter  than  that  of  the  uracil* diaminopurine 
duplex  that  has  three  hydrogen  bonds  (UF*D,  Kb  =  820  nM). 
Similar  tight  binding  was  observed  for  the  duplex  with  an 
abasic  site  opposite  UF  (UF*3>,  £d  =  2.4  nM),  providing 
additional  evidence  that  complete  ablation  of  hydrogen 
bonding  leads  to  a  significantly  increased  binding  affinity.3 


3  The  observed  KD  for  the  UF*d>  duplex  (2.5  nM)  predominantly 
reflects  binding  of  UDG  to  UF,  and  not  to  <£,  because  the  KD  for  abasic 
DNA  is  much  weaker  (70  nM  for  T*d>  and  100  nM  for  4>*A). 
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Table  1:  Thermodynamic  Parameters  for  UDG  Binding,  DNA  Melting,  and  Permanganate  Accessibility  of  Destabilized  Base  Pairs 


duplex  Ky  (nM) _ AGbind  (kcal/mol) _ A//meit  (kcal/mol) _ A5mcit  (cal  mol  1  K  *) _ AGmeit°  (kcal/mol) _ log  & 


UF-D 

820  ±  90 

-8.3  ±0.1 

58.9  ±  1.8 

UF-A 

300  ±50 

-8.9  ±0.1 

39.6  ±  4.0 

UF-N 

77  ±11 

-9.7  ±0.1 

40.5  ±2.0 

UF-G 

57  ±  6 

-9.9  ±0.1 

34.0  ±  2.4 

UF-4>3 

2.5  ±0.4 

—  1 1 .7  ±0.1 

4.3  ±0.1 

UF-M 

0.24  ±  0.03 

— 13.1  ±  0.1 

NDC 

°  Calculated  at  298  K.  b  S  is  the  relative  sensitivity  of  a  T*X  base  pair 


0.16  ±0.01 

10.9  ±2.3 

-1.4  ±0.14 

0.11  ±0.01 

7.1  ±5.2 

-1.0  ±0.06 

0.10  ±0.01 

9.2  ±  2.5 

ND* 

0.10  ±0.01 

5.6  ±  3.2 

-0.80  ±0.03 

0.01  ±0.01 

0.3  ±  0.3 

-0.22  ±0.01 

NDC 

NDC 

-0.36  ±0.01 

to  oxidation  by  KMnCU  (see  the  legend  of  Figure  6). c  Not  determined. 


These  results  indicate  that  removal  of  three  hydrogen  bonds 
can  enhance  specific  recognition  by  up  to  4.8  kcal/mol. 

In  large  part,  these  binding  effects  reflect  the  thermo¬ 
dynamic  properties  of  the  free  damaged  site  because  these 
discrete  base  pair  perturbations  are  not  expected  to  affect 
interactions  between  UDG  and  the  DNA.4  This  conclusion 
is  supported  by  inspection  of  the  crystal  structures  of  UDG 
complexed  with  substrate  analogues,  which  show  that  UDG 
does  not  make  any  interactions  with  the  base  that  opposes 
the  damage  site,  or  with  the  undamaged  strand  (36).  Thus, 
recognition  solely  involves  the  extrahelical  deoxyuridine  and 
not  other  specific  features  of  the  base  pair  or  duplex.  To 
further  establish  that  all  of  the  UF*X  duplexes  used  here  attain 
the  same  bound  state,  and  that  the  observed  effects  on  binding 
largely  arise  from  the  properties  of  the  free  DNA,  we 
measured  the  tryptophan  fluorescence  quenching  upon  bind¬ 
ing  of  each  duplex  (Figure  4B).  Previous  work  has  shown 
that  the  quenching  of  UDG  tryptophan  fluorescence  upon 
specific  DNA  binding  is  a  sensitive  measure  of  an  induced 
fit  conformational  change  in  UDG  that  is  required  to  achieve 
the  final  productive  conformation  with  a  flipped-out  uracil 
(5).  Within  the  errors  of  these  measurements,  all  of  the 
duplexes  produced  the  same  magnitude  of  fluorescence 
quenching,  indicating  that  the  same  bound  conformations 
were  attained  for  all. 

Thermodynamic  Stabilities  of  IF'X  Duplexes.  We  then 
determined  the  energetic  effects  of  this  series  of  site-specific 
base  pair  disruptions  on  the  thermodynamic  parameters  for 
DNA  melting  using  differential  scanning  calorimetry  (DSC, 
Figure  5).  With  DSC,  one  can  measure  the  enthalpy  (AHmc\t) 
and  entropy  (ASmeit)  of  the  melting  transition  directly,  and 
unlike  optical  methods,  it  is  insensitive  to  the  mechanism  of 
duplex  melting  (37).  The  free  energy  of  duplex  dissociation 
(AGmeit)  at  any  temperature  can  then  be  simply  calculated 
from  AZ/meit  and  A5meit,  using  the  relationship  AGmeit  = 
AZ/ineit  “  TASmc\t ,  because  the  heat  capacity  of  the  duplex 
and  that  of  the  single  strands  are  equal  (A Cp  &  0)  (38,  39). 
As  expected,  UF-X  duplexes  exhibited  decreasing  transition 
enthalpies  in  the  A//meit  range  of  58.9—4.3  kcal/mol  as  the 
number  of  hydrogen  bonds  was  decreased.  The  complete 
thermodynamic  parameters  for  duplex  melting  are  reported 
in  Table  1. 

Correlation  of  UDG  Binding  Affinity  with  Damaged  Base 
Pair  Stability.  To  quantitatively  evaluate  the  impact  of 
damaged  base  pair  disruption  on  UDG  affinity,  the  differ- 


4  The  observation  of  linear  free  energy  relationships  between  the 
thermodynamic  parameters  of  the  free  duplex  DNA  and  the  overall 
free  energy  of  DNA  binding  is  not  negated  if  different  energetic 
interactions  exist  between  the  enzyme  and  each  bound  duplex.  However, 
the  observed  slopes  will  reflect  the  relative  effects  of  perturbing  the 
duplex  in  the  free  and  bound  state. 


T(°C) 

Figure  5:  Differential  scanning  calorimetry  (DSC)  of  the  UF*D 
duplex.  A  20  pM  DNA  duplex  solution  in  10  mM  sodium  phosphate 
(pH  7.5)  and  60  mM  NaCI  was  placed  inside  of  the  DSC  sample 
cell,  and  the  change  in  the  heat  capacity  of  the  solution  (AC,,)  was 
monitored  as  it  was  warmed  from  20  to  95  °C  at  a  rate  of  60  °C/h. 
The  thermodynamic  parameters  of  duplex  dissociation  were  ex¬ 
tracted  as  described  in  Experimental  Procedures. 

ences  in  binding  free  energies  (AAGbind,  Figure  1)  were 
plotted  against  the  changes  in  transition  enthalpies  (AA/fmeit) 
and  entropies  (-7AASmeit)  (Figure  6A,B).  In  this  analysis, 
the  difference  energies  are  relative  to  the  duplex  with  three 
hydrogen  bonds  (UF*D),  and  the  value  of  -TAASme]t  was 
calculated  at  298  K.  AAGbind  was  found  to  increase  linearly 
as  the  transition  enthalpy  and  entropy  decrease:  AAHme\t 
(slope  =  0.064  ±  0.01,  r2  =  0.937)  and  —TAASme\t  (slope 
=  0.080  ±  0.01 1,  t2  —  0.949).  As  required  from  these  linear 
correlations,  a  strong  correlation  with  AAGmeit  (slope  =  0.3 
±  0.1,  r2  =  0.839)  was  also  observed  (data  not  shown).  The 
implications  of  these  correlations  are  discussed  below. 

Correlation  of  UDG  Binding  Affinity  with  Damaged  Base 
Accessibility.  The  above  thermodynamic  correlations  suggest 
that  these  destabilized  duplexes  might  exhibit  an  increase  in 
the  number  of  dynamic  fluctuations  that  promote  extrahelical 
states  of  UF  at  temperatures  well  below  the  duplex  melting 
temperature.  Since  the  dynamic  accessibility  of  the  damaged 
base  is  another  factor  that  could  enhance  its  recognition  by 
DNA  glycosylases,  it  was  of  interest  to  measure  the  relative 
accessibility  of  each  destabilized  base  pair,  and  correlate  this 
parameter  with  UDG  binding  affinity. 

To  explore  this  question,  a  potassium  permanganate 
(KMnC>4)  sensitivity  assay  was  employed  (32, 35,  40).  Since 
sites  of  pyrimidine  oxidation  are  susceptible  to  strand 
cleavage  under  basic  conditions,  they  can  be  detected  as 
fragments  using  polyacrylamide  gel  electrophoresis  (Figure 
7A)  (41).  In  these  studies,  UF  was  first  replaced  with  thymine 
(T)  because  of  its  poor  oxidation  reactivity  arising  from  its 
electron  deficient  5,6-double  bond  (data  not  shown).  This  is 
a  very  conservative  change,  as  a  T*X  base  pair  will  have 
base  pairing  strength  and  geometry  nearly  identical  to  those 
of  a  UF*X  base  pair  (42,  43).  As  shown  in  Figure  7 A, 
decreasing  the  base  pair  strength  leads  to  an  increase  in  the 
sensitivity  of  the  T*X  base  pair  to  permanganate  oxidation, 
and  also  the  invariant  T- A  base  pair  three  nucleotides  away, 
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Figure  6:  Correlations  between  UDG  binding  affinity  for  each 
UF*X„  duplex  (AAGbind)  and  the  thermodynamic  parameters  for 
duplex  melting.  The  numbers  of  hydrogen  bonds  in  the  UF*X„  base 
pair  are  indicated.  (A)  AAGbind  vs  AA//meit  (slope  =  0.064,  r2  = 
0.937).  Only  an  upper  limit  for  the  enthalpy  of  melting  for  the  UF* 
M  duplex  was  obtained  (arrow).  (B)  AAGbind  vs  —TAAS^u  (slope 
=  0.080,  r2  =  0.949). 

indicating  that  even  very  conservative  changes  to  the  T\X 
base  pair  can  influence  the  dynamics  of  neighboring  base 
pairs  in  the  duplex  (44).  A  plot  of  log  KD  against  log(relative 
KMn04  sensitivity)  shows  a  linear  correlation  (Figure  7B), 
establishing  that  increasing  base  accessibility  at  temperatures 
well  below  the  Tm  value  has  a  strong  positive  effect  on 
binding  affinity.  As  an  important  control,  the  amount  of 
oxidized  product  was  found  to  increase  linearly  with  respect 
to  time  and  concentration  of  KMn04  (data  not  shown). 
Therefore,  the  differences  in  the  sensitivity  of  these  duplexes 
to  oxidation  directly  reflect  the  unfavorable  dynamic  pre¬ 
equilibrium  for  exposure  of  the  thymidine  base  prior  to 
reaction  with  KMn04. 

DISCUSSION 

Thermodynamic  Framework  for  Active  and  Passive  Base 
Flipping.  While  it  is  clear  from  structural  and  spectroscopic 
studies  that  DNA  glycosylases  bind  their  cognate  damaged 
base  in  an  extrahelical  conformation  (5,  11)>  the  pathway 
by  which  the  damaged  base  is  flipped  out  of  the  DNA  duplex 
and  placed  inside  of  the  enzyme  active  site  remains  poorly 
defined.  In  one  model,  DNA  glycosylases  passively  capture 
damaged  bases  that  are  transiently  extrahelical  (Figure  8) 
(14,  45).  According  to  this  view,  the  DNA  glycosylase  does 
not  lower  the  activation  energy  or  equilibrium  for  damaged 
base  flipping,  but  instead  relies  on  the  increased  extrahelical 
propensity  of  damaged  bases  to  enhance  bimolecular  en- 
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Figure  7:  (A)  Sensitivity  of  thymine  in  single-stranded  and  duplex 
DNA  to  oxidation  by  potassium  permanganate.  Samples  of  5'-32P- 
labeled  single-stranded  or  T*X„  duplex  DNAs  were  reacted  with 
2.5  mM  KMn04  ([O])  for  3  min.  After  the  oxidized  strands  had 
been  cleaved  with  piperidine,  the  reaction  mixtures  were  run  a  19% 
denaturing  polyacrylamide  gel,  and  the  radioactivity  of  each  band 
was  quantified  with  a  phosphorimager.  The  position  on  the  gel  of 
the  full-length  DNA  is  marked  with  a  black  wedge.  (B)  Correlation 
between  UDG  binding  affinity  and  KMn04  sensitivity  (slope  = 
—2.37  dt  0.270,  r2  =  0.963).  The  relative  sensitivity  of  a  thymine 
in  a  T*X„  base  pair  is  defined  as  [(T18  -  7»*gd)/(/totai  -  ^duplex  x 

lOO]/^8  -  /*gd)/(/total  _  ^bkgdyingle— stranded  x  100],  where  F*  is 

the  intensity  of  the  band  corresponding  to  oxidation  of  the  T  in  the 
T*X„  base  pair,  7bkgd  is  the  background  correction,  and  Z01®1  is  the 
sum  of  all  of  the  intensities  of  the  bands  in  a  given  lane.  The 
numbers  of  hydrogen  bonds  in  the  T*X„  duplex  are  indicated. 
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Figure  8:  Thermodynamic  model  showing  the  energetic  equiva¬ 
lence  of  the  passive  and  active  base  flipping  mechanisms  (see  the 
text).  The  equivalence  of  the  pathways  is  illustrated  by  the 
thermodynamic  box,  which  requires  that  l/XDapp  =  XjnbindXenzfijP  = 
XnonflipXoutbind. 


counter.  In  the  alternative  view,  DNA  glycosylases  actively 
flip  out  their  cognate  lesions  by  destabilizing  the  damaged 
site  in  an  initial  encounter  complex.  Active  base  flipping 
may  occur  by  stabilization  of  high-energy  intermediate 
conformations  on  the  base  flipping  pathway,  or  by  the  use 
of  mechanical  forces  to  propel  the  base  from  the  duplex 
(Figure  8)  (5,  46).  It  is  important  to  point  out  that  passive 
and  active  flipping  pathways  cannot  be  distinguished  by 
thermodynamic  measurements  alone .  This  conclusion  is 
required  because  =  XnonfliP^outbind,  as  shown  in 
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Figure  8.  Only  an  assessment  of  the  kinetic  competence  of 
each  pathway  can  reveal  whether  passive  or  active  base 
flipping  is  the  major  route  taken  for  a  given  system. 

Specific  Recognition  and  Conformational  Freedom  of 
Damaged  Sites.  We  have  found  that  the  affinity  of  UDG  for 
a  specific  site  is  linearly  dependent  on  the  enthalpy  of  duplex 
dissociation  (Figure  6A).  The  linear  relationship  between 
UDG  affinity  and  A//meit  confirms  the  notion  that  weakened 
base  pairing  decreases  the  enthalpy  of  duplex  melting  and, 
consequently,  increases  UDG  binding  affinity  (Figure  6 A). 
However,  this  finding  alone  provides  little  insight  into  the 
physical  mechanism  by  which  decreases  in  base  pair  enthalpy 
lead  to  enhanced  binding.  The  small  slope  of  the  correlation 
suggests  that  the  higher  affinity  of  UDG  for  destabilized 
damaged  sites  does  not  simply  arise  from  the  reduced 
enthalpic  cost  of  breaking  a  destabilized  base  pair  during 
base  flipping,  and  that  other  energetic  influences  must  be  at 
work.4  It  should  be  stated  that  the  experiments  presented  here 
are  explicitly  designed  to  probe  the  enthalpic  and  entropic 
contributions  of  the  damaged  site  alone  to  specific  recogni¬ 
tion.  The  total  enthalpy  and  entropy  of  binding  (including 
the  enzyme,  DNA,  and  solvent)  are  not  evaluated  in  any  of 
the  current  experiments,  and  in  fact,  these  measurements  are 
not  required  for  the  conclusions  presented  below. 

The  plot  of  AAGbind  versus  -TAASmc\t  reveals  that  there 
is  an  equally  significant  correlation  between  changes  in 
duplex  entropy  and  binding  affinity  (Figure  6B).  The 
parameter  “7AASmeit  likely  reflects  the  increased  confor¬ 
mational  flexibility  of  the  destabilized  base  pairs  in  the 
duplex  DNA,  because  the  entropy  differences  of  the  dis¬ 
sociated  single-stranded  DNAs  in  the  melting  experiments 
should  be  similar,  given  the  conservative  changes  in  these 
substrates.  If  we  assume  this  physical  interpretation  for  the 
entropy  changes  between  these  DNA  constructs,  the  cor¬ 
relation  suggests  that  increased  flexibility  of  the  base  pair 
produces  conformational  states  that  are  productive  for  UDG 
binding.  The  conclusion  that  extrahelical  conformational 
states  are  produced  is  supported  by  the  correlation  between 
UDG  binding  affinity  and  KMn04  sensitivity  (Figure  7B), 
which  reflects  the  dynamic  equilibrium  of  the  base  between 
an  inaccessible  and  permanganate  accessible  state  (i.e.,  an 
extrahelical  exposed  conformation).  A  reasonable  interpreta¬ 
tion  of  these  combined  findings  is  that  enthalpic  destabiliza¬ 
tion  of  the  base  pair  allows  increased  conformational 
flexibility,  producing  extrahelical  conformers,  some  of  which 
favor  enzyme  binding. 

The  effects  of  base  pair  enthalpy  and  entropy  on  the 
extrahelical  conformational  distributions  that  may  promote 
base  flipping  are  depicted  in  Figure  9A— C.  In  these  panels, 
the  probability  of  an  extrahelical  conformation  is  plotted 
against  the  backbone  pseudodihedral  angle  (\p)  of  the 
deoxyuridine  nucleotide,  defined  as  indicated  in  Figure  9. 
The  angle  xp  has  been  previously  used  in  computational 
studies  to  describe  the  pathway  for  base  flipping  (47),  and 
is  used  here  because  of  its  simple  representation  of  the  base 
flipping  trajectory,  although  none  of  the  arguments  depend 
on  this  formalism.  Using  this  nomenclature,  a  xp  of  10° 
reflects  the  fully  base  paired  state  and  a  xp  of  180°  reflects 
the  fully  extrahelical  state.  In  free  DNA  (Figure  9A),  the  U 
nucleotide  in  the  stable  base  pair  with  three  hydrogen  bonds 
(U*D)  should  be  tightly  centered  around  an  average  confor¬ 
mation  with  a  xp  of  10°  (red  curve),  while  the  unstable 
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Figure  9:  Enthalpic  destabilization  of  the  U*X„  base  pair  leads  to 
extrahelical  conformations  that  promote  binding.  The  probability 
of  an  extrahelical  conformer  as  a  function  of  the  pseudodihedral 
angle  xp,  which  is  a  measure  of  the  progress  along  the  base  flipping 
reaction  coordinate,  is  indicated  (see  the  text  and  ref  47).  Using 
this  nomenclature,  a  fully  stacked  base  pair  has  a  xp  value  of  M0°, 
whereas  the  fully  extrahelical  conformation  has  a  xp  value  of  1 80°. 

(A)  Enthalpic  destabilization  of  the  base  pair  in  the  free  DNA  leads 
to  extrahelical  conformations  that  promote  passive  base  flipping. 

(B)  Base  pair  destabilization  can  also  affect  an  active  mechanism 
in  which  UDG  forms  an  initial  encounter  complex  with  the  DNA 
in  which  the  base  is  not  yet  fully  extrahelical  (5,  46,  48).  In  this 
case,  enzyme  binding  energy  is  used  to  destabilize  the  U*D  base 
pair,  allowing  it  to  achieve  extrahelical  states  already  available  to 
the  U*<I>  duplex  due  to  its  intrinsic  instability.  (C)  A  hypothetical 
distribution  of  extrahelical  conformers  for  the  final  UDG-DNA 
complex.  The  enzyme  has  fully  stabilized  the  flipped-out  uracil, 
and  the  distribution  of  conformational  states  for  the  U*D  and  U*<b 
duplexes  is  narrowly  focused  around  a  xp  of  180°  (36,  47). 


construct  with  no  hydrogen  bonds  (U*3>)  should  be  broadly 
centered  around  a  xp  value  further  along  the  flipping  reaction 
coordinate  (green  curve).  Thus,  enthalpic  destabilization  of 
the  base  pair  in  a  passive  mechanism  leads  to  conformations 
that  facilitate  binding  of  the  enzyme  (see  Figure  8,  counter¬ 
clockwise  pathway). 

For  the  alternative  active  base  flipping  mechanism  (Figure 
9B),  binding  of  UDG  in  an  initial  encounter  complex  can 


4194  Biochemistry,  Vol  43,  No.  14,  2004 

alter  the  DNA  structure  such  that  the  more  stable  U*D  base 
pair  is  also  destabilized  (Figure  9B,  red  curve)  (5,  46,  48). 
In  this  initial  destabilized  complex,  the  partially  extrahelical 
uracil  may  assume  a  similar  average  conformation,  and  broad 
distribution,  for  both  the  U-D  and  U •<$  constructs.  However, 
flipping  of  U  in  the  U*D  construct  requires  a  greater  amount 
of  binding  energy  to  overcome  the  enthalpic  barrier  to  base 
flipping,  resulting  in  weaker  binding  of  the  U*D  duplex 
compared  to  that  of  the  duplex.  Finally,  in  the  Michaelis 
complex,  in  which  the  base  is  frilly  extrahelical  (Figure  9C), 
both  the  U'D  and  U-d>  duplexes  assume  the  same  average 
conformation  and  small  conformational  distribution  which 
are  enforced  by  the  strong  interactions  between  the  enzyme 
and  uracil  (36,  46).  Although  this  mechanism  implies  a 
reduction  in  the  conformational  flexibility  of  destabilized 
base  pairs  such  as  the  U*0  duplex  upon  formation  of  the 
Michaelis  complex,  this  expected  unfavorable  entropic 
contribution  to  the  overall  free  energy  of  binding  may  be 
paid  for  by  the  even  larger  enthalpic  benefit  of  base  pair 
destabilization  for  such  conformationally  flexible  substrates 
(Table  1). 

CONCLUSIONS 

We  have  shown  that  the  maximum  gain  in  specific 
recognition  by  UDG  arising  from  destabilization  of  the 
damaged  base  pair  is  4.8  kcal/mol,  and  that  this  energetic 
effect  likely  arises  from  increasing  the  population  of  extra¬ 
helical  states  that  promote  binding.  Since  the  catalytic 
specificity  of  UDG  for  uracil  as  opposed  to  other  normal 
bases  has  been  estimated  to  be  at  least  8.3  kcal/mol  in  vitro 
(5),  the  energetic  contribution  of  destabilized  damaged  sites 
to  specific  ground  state  binding  can  be  significant,  at  least 
in  this  model  system.  However,  the  in  vivo  substrates  of 
UDG  consist  of  U*A  or  U*G  base  pairs,  which  are  not 
significantly  destabilized  compared  to  other  normal  base 
pairs.  Thus,  spontaneous  base  flipping  at  damaged  sites  is 
not  a  viable  mechanism  for  accounting  for  the  specificity  of 
UDG  in  vivo.  The  remaining  specificity  of  UDG  must  be 
attributed  to  strong  transition  state  interactions  that  can  be 
induced  only  by  actively  flipping  the  uracil  base  into  the 
active  site  (5,  49).  Although  unimportant  for  UDG,  damaged 
site  instability  could  contribute  significantly  to  specific 
recognition  by  repair  enzymes  that  act  on  intrinsically 
unstable  base  pairs  such  as  C^-methylguanine  (45),  hypox- 
anthine  (50,  51),  A^-methyladenine,  and  A^-methylcytosine 
(52,  53). 
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